Distillation of Multipartite Entanglement by Complementary Stabilizer Measurements 
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We propose a scheme of multipartite entanglement distillation driven by a complementary pair 
of stabilizer measurements, to distill directly a wider range of states beyond the stabilizer code 
states (such as the Greenberger-Horne-Zeilinger states). We make our idea explicit by constructing 
a recurrence protocol for the 3-qubit W state -^=(|001) + 010) + 100) ) . Noisy W states resulting 
from typical decoherence can be directly purified in a few steps, if their initial fidelity is larger than 
a threshold. For general input mixed states, we observe distillations to hierarchical fixed points, i.e., 
not only to the W state but also to the 2-qubit Bell pair, depending on their initial entanglement. 

PACS numbers: 03.67.Mn, 03.67.Pp, 03.65.Ud, 42.50.-p 



Introduction.- We recognize, through a decade of re- 
search, that entanglement is indispensable to execute 
quantum information processing (QIP), such as quantum 
computation and multi-party quantum communication. 
A persistent challenge is to maintain multipartite entan- 
glement against decoherence. In this Letter, we enlarge 
the present applicability of a key technique, entangle- 
ment distillation 0> Q] , to genuine multipartite entangled 
states called the W states 

The W state, -^(|0 . . . 01) + 10 . . . 10) + • • • + 110 . . . 0)), 
i.e., the equal superposition of "single-excitation" basis 
vectors in n qubits is a tolerant resource against decoher- 
ence and loss of qubits. It is quite robust 0,O|, because 
it can be compared to a symmetric web consisting of 
only pairwise entanglement @, @. This state is also a 
Dicke state of the total spin operators J 2 and J z with 
eigenvalues f (^ + 1) and § — 1 respectively, due to the 
permutation symmetry. Thus the W state is often avail- 
able more easily than the Greenberger-Horne-Zeilinger 
(GHZ) state. The 3-qubit W state has been created in 
optical systems Q and ion traps Q , and can be prepared 
according to several proposals in coupled quantum dots, 
critical spin chains, etc. Furthermore, the W state can 
be said to be essentially different from most multipartite 
entangled states known in the applications to QIP, in the 
following sense. 

Basic software techniques to circumvent decoherence 
have been proposed, and implemented experimentally for 
prototypes. For example, entanglement distillation (or, 
purification) 0,0 is a tool to extract high-fidelity entan- 
gled states from a larger ensemble of noisy ones. Quan- 
tum error correction codes 0, flfl ] are a way to protect 
entanglement from small numbers of errors. The latter 
can be formulated in terms of the stabilizer, i.e., as si- 
multaneous eigenspaces of commuting "multifocal" Pauli 
operators. Note that, if all stabilized eigenspaces are 1- 
dimensional state vectors, they are often called stabilizer 
states or graph states (up to local unitaries) 0, 0] . 
The Calderbank-Shor-Steane (CSS) code Q is defined by 



the stabilizer group which consists of only two kinds of 
generators: multifocal bit-flip X operators and phase-flip 
Z operators. In fact, beyond the "bipartite" distillation 
protocols for the Bell pairs, direct distillation of 

multipartite entanglement is so far possible just for the 
CSS stabilizer (or, two-colorable graph) states by the pro- 
tocol in Ref s. Il4lllffl . which extended earlier results for 
GHZ states [l|l [l7| . Since the W state is not a stabilizer 
state, there has been no protocol to distill it directly. 

Main idea.- We propose an entanglement distilla- 
tion protocol that extracts directly a multipartite non- 
stabilizer (non-graph) state, specifically the 3-qubit W 
state. Our idea is to apply local measurements of the 
stabilizer (whose nonlocal counterpart acting at differ- 
ent parties stabilizes the target state), assuming that the 
target state belongs to a basis of equivalent entangled 
states. Note that such a basis, similar to the Bell basis, 
exists for a wider range of multipartite states than stabi- 
lizer states. We need n copies of the input state for the 
n-qubit case, to apply stabilizer measurements locally. 
In this manner, we can improve the fidelity, and attain 
the target state as a fixed point of the protocol. Note 
that if the target state is not the stabilizer state, local 
depolarization or twirling (over the single copy) which 
keeps the target state invariant seems impossible in gen- 
eral. Thus, we do not make the mixed states diagonal, 
i.e., a classical mixture of the basis states. It implies 
we cannot reduce the task to a "classical problem" that 
consists in extracting entropy from the binary strings of 
the stabilizer eigenvalues, as is possible by bilateral CNOT 
operations in all the known protocols. Nevertheless, by 
virtue of complementary stabilizer measurements which 
exchange the amplified components, our protocol works 
without local depolarization. The feature is favorable in 
efficiency, and analogous to the Oxford protocol 

Direct distillation of multipartite entanglement has 
several potential merits. In the case of CSS states such as 
the GHZ states, multipartite distillation was shown to be 
more efficient than the bipartite strategy which consists 
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of distillation of Bell pairs and their connection [Til . 
Under imperfect operations, the achievable fidelity can 
be higher 

003. Also, the threshold for distillability 
may be tighter than that by indirect methods. 

W basis and its stabilizer group.- To make our idea 
explicit, we construct a recurrence protocol for the 3- 
qubit W state \W om ) = ^(|001) + |010) + |100)) ABC , 
distributed over Alice, Bob, and Carol. We denote the 
Pauli matrices, operating on the j-th qubit at the party 
I, as Aj, Yj, and Z l - along with the identity t^. To 
distinguish the non-local tensor structure of the multiple 
Hilbert spaces controlled by different parties and the local 
tensor structure at a single party, we use the superscripts 
I as the non-local indices, and the subscripts j as the 
local ones. We define the Hadamard operation by H = 
-^(X+Z), the 2-qubit swap operation by swap : \kk') 1— » 

\k'k) (k,k' = 0,1) in the computational basis, and a 3- 
qubit unitary operation V, which leaves 1 000) and 1 1 1 1 ) 
unchanged, but exchanges the others in such a way that 
|001) <-> 1 110) , |010) <-> 1 101) , and 1 100) <-► |011). 

Let us introduce a complete orthonormal basis, called 
the W basis here, where each basis state \W klk2k3 ) 
(k x ,k 2 ,k 3 — 0,1) has entanglement equivalent to the 
W state |W 000 ). This is because basis states transform 
into each other by the local unitary operations in Ta- 
ble [I] The W basis can be obtained from the compu- 
tational basis acted on by the 3-qubit unitary operation 
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\ w k t k 2 k^ = u wh&sis \kik 2 k 3 ). It is convenient to identify 
the stabilizer for the W basis. To satisfy K 3 \W klk2k3 ) — 
(— \) kj ' \W klk2ks ), three generators Kj are determined as 

K (ABCL i {2 X A X B Z c + 2Y A Z B Y C + Z A t B l c ), 
K (ABC L i {2Z a x b x c + 2Y A Y B Z C + t A Z B t c ), (1) 
K (ABCL i { 2x A Z B X c + 2Z A Y B Y C + 1 A 1 B Z C ). 

We emphasize, by the superscript (ABC), that the sta- 
bilizers are not local. Note that later we measure locally 
the stabilizers, which will be denoted, e.g. for Alice, as 
K[ A) = \{2X A X A Z A + 2Y A Z A Y A + Z A 1 A 1 A ), and all 
3 qubits specified by the subscripts belong to Alice. The 
stabilizer group consists of eight commuting elements, 
{1, K X ,K 2 , K 3 , K X K 2 , K X K 3 , K 2 K 3 , K X K 2 K 3 }, where 



K X K 2 = \(2\ X X 2 X 3 + 2Y X \ 2 Y 3 - Z x Z 2 t 3 
K X K 3 = \(2X x X 2 t 3 + 2t x Y 2 Y 3 - Z x t 2 Z 3 
K 2 K 3 = \(2X x t 2 X 3 + 2Y x Y 2 t 3 - t x Z 2 Z 3 



(2) 



K X K 2 K 3 



-Z X Z 2 Z 3 



W state distillation protocol. - Our protocol consists of 
two subprotocols: £P and its dual 3? . In both, three in- 
put copies are mapped into one output copy to define a 
simple recurrence. Generally, any mapping to a smaller 
subsystem can be considered. We assume, without loss of 
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TABLE I: The 3-qubit W basis. Local unitary operations that 
map |W 000 ) to \ W klk2k3 } are shown in the third column. 



generality, that the |W 000 )(W 000 | component of the input 
mixed states p is the largest among the diagonal elements 
in the W basis (otherwise we can relabel the computa- 
tional basis by a local unitary operation in Table [IJ. We 
define the fidelity F of p by F = (W OOQ \p\W om ). 

Protocol 1) Every party (/ = A, B, or C) applies 
the local measurement of two stabilizers K^p and 



K X K% over the input state 7 (= p 
and obtains the 2-bit outcomes = [m x ' , nia']- 2) In 
forming their outcomes by two-way classical communica- 
tion, parties select coincident outcomes = = 
m( c ) = [0,1] (= 1), [1,0] (= 2), or [1,1] (= 3). Oth- 
erwise they discard three copies. 3) For the coincident 
outcomes, each party transforms locally her/his state into 
a 1-qubit subsystem by the following "majority rule" . 
If m« = 1, Pi : I Wboi)' >-> \0) 1 , \W XW ) L ^ if 
mW =2, Pi: \W 0W ) 1 h- \0) 1 , \W WX ) 1 h- |1)'; and if 
m« = 3, Pi : |VL 10 o)' ^ W, \W mx ) 1 h-> |1)'. 

Mathematically, the stabilizer measurement M (,% 
the party I is written by, 
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with the completeness condition 

Note that M ( acts on 3 qubits of the party I, and it 
is a projector to the local W basis vectors, for example if 
m« = 1, MjP = \W 00X ) 1 (W001 1 + |Wiio)'(Wno|. By the 
selection of desired coincident outcomes m, & maps the 
input state 7 (= p^ 3 ) to the one-copy state p' given by 

p' = J2 PM^PM^PM£hPM^PM^PM^\ 

m=l,2,3 

(4) 

with the success probability tr(p'). We normalize the 
state as p out = p'/tr(p') for the next recurrence step. 

Before describing the whole protocol including we 
illustrate analytically how & works. Suppose the perfect 
W state is distributed to three parties, but suffers typical 
decoherence as described by the local dephasing channel 
V l (p) = \{{l+p)p+(l-p)Z l pZ l ) with the same channel 
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initial fidelity 



FIG. 1: The distillation curve (left) of 0* for locally dephased 
W states, and the yield (right) after achieves F > 0.99. 
Note that the region of F G [0, ^) is not physical. 

reliability p, G [0, 1]. Three parties initially share a noisy 
W state a(F) = V A V B V c {\W om ){W mQ \), which is not 
diagonal in the W basis, but is parametrized uniquely by 
F = |(1 + 2p 2 ) G [g,l]< A straightforward calculation 
shows that 3? maps three copies tr(F)® 3 to one copy 
a(F') with the higher fidelity F' such that 

p l= f^ + ^F(l-Ff + ^{l-Ff 

- || F 3 + 1 F 2 (1 _ F)+ 2_ F{1 _ F) 2 + (1 _ p )3 ■ 

(5) 

Eq. (J5J suggests a recurrence seen in the distillation curve 
of Fig. n Since F lies in [g, 1], we prove analytically 
that F = 1, corresponding to the W state, is the attrac- 
tive fixed point and F = | is the repulsive one. We find 
that for any locally dephased W state except F = h, 3? 
restores it with a few steps. Indeed, this threshold coin- 
cides with a necessary condition for distillability by the 
partial transposition criterion 19, 20] . Since the mixed 
state T l (a(F)), partially transposed for any party I (i.e., 
bipartition), has a negative eigenvalue only for F > i, 
there is no chance to distill entanglement in F — i. In 
Fig. n the yield (i.e., the ratio of the number of surviv- 
ing copies to that of used copies) of after F reaches at 
least 0.99 is also shown. The "stairs" of the yield come 
from the difference in the number of recurrence steps. 

For more general noises, we need 3^ which has the 
similar structure as 3^ but employs complementary ob- 
servables Rf = k^Kfk 1 , where A' = H[H^H 1 3 swap[ 3 . 
Two measurement bases \Wk 1 k 2 k 3 ) in 3* and \Wk'k' 2 k'.) — 
h'\Wk' h^k's) m & are complementary (also called mutu- 
ally unbiased [111), i.e., \{W klk2k3 \W k > ik > 2 k> 3 )\ 2 = |. 

Dual Protocol 3?: 0) Every party (/ = A, B, or C) 
applies V 1 to change the local computational basis. The 
input 7 is modified to 7 = V A V B V C pf*V A W B ^V c ^ . 1) 
She/he applies the local measurement of two dual stabi- 
lizers K2 and K 3 on 7, and obtains the 2-bit 
outcomes m!'l 2) By two-way classical communication, 
parties select the coincident outcome = m( B ) = 

m( c ) = [0,0] (= 0). 3) Each party transforms locally 
the state into a 1-qubit subsystem in the manner oppo- 



locally dephased W states locally depolarized W states 




FIG. 2: Noisy W states subjected by the local dephasing (left) 
or local depolarizing (right) channel can be retrieved by SP 
and 0, if F is initially larger than | or 0.48, respectively. 
Note that this is actually accomplished by alone for the 
local dephasing case. 

site to 9\ P l Q : IIFqoo}' H\l) 1 , \W ln ) 1 ~ H\0} 1 . 

In brief, in 3^, we replace all operators in Eqs. © and 
QJ by their "barred" dual operators. The complete dis- 
tillation procedure for general mixed states consists of the 
sequential application of either 3^ or 3^ where, in every 
recurrence step, we select one of the subprotocols which 
gives the higher fidelity in the output. There seems to 
be no simple formula for the sequence, but we can de- 
termine the sequence if we know the initial input state p 
before distillation, e.g., by state tomography. Also note 
that although this combination of 3? and & can reach 
numerically the region where F > 0.999 , precisely speak- 
ing F = 1 is the fixed point of 3? but is not that of 0*. 

Hereafter, we show that, under the sequential ap- 
plication of 3? and the W state can be distilled 
from arbitrary mixed states if, roughly speaking, F is 
sufficiently large. First, consider another typical deco- 
herence such as the local depolarizing channel (white 
noise) £ l (p) = W + ^{p + X 1 pX l + Y l pY l + Z l pZ l ), 
and the input state p in (F) = £ A £ B £ C (\W mo ){W 000 \) 
with F = i(3 + p + 9p 2 + 11/i 3 ) G [|, 1]. Although 
the locally depolarized W state does not remain in the 
same form under our protocol, we can still determine a 
threshold for distillability. As seen in Fig. |2 if initially 
F > 0.48, we distill the W state, and otherwise we have 
an undistillable mixed state x = ji^)^ + l^'K^'Di 
where \ip) = |(|001) + |010) + 100) - | 111)) and \ip') = 
5HOOO) + |011) + 1 101) + 1 110)), as another fixed point 
with F = |. This threshold is stricter than the neces- 

o 

sary condition F > 0.36 by the partial transpose criterion 
|l9l l20j . Note that the progress of the protocol is not de- 
scribed by a single parameter, and F is not monotonic 
any more. A nonmonotonic behavior of F was also seen 
in the bipartite distillation without depolarization |l3j . 
However, as a long-term behavior, F is increasing for the 
distillable cases and can be used for visualization of the 
progress. 

Next, we consider randomly generated input mixed 
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random states (F ~ 0.70) random states (F ~ 0.50) 




FIG. 3: The average fidelity and its standard deviation fol- 
lowed toward each fixed point, for (in total) 10,000 randomly 
generated initial mixed states with F £ 0.70 ± 0.01 (left) or 
0.50 ±0.01 (right). 



states (under the Hilbert-Schmidt measure 22J), and will 
observe numerically hierarchical distillations not only to 
the 3-qubit W state, but also to a 2-qubit Bell pair. This 
is surprising, since it implies that we can distill a non- 
stabilizer state and a stabilizer state by the same proto- 
col. In Fig. |3 for 10,000 random mixed states with the 
initial fidelity F fixed close to 0.70 or 0.50, we display 
the average fidelity and its standard deviation for each 
set of samples reaching the same fixed point. When F is 
sufficiently large, such as F ~ 0.70, the branch to the W 
state is dominant. More than 99 percent of the states fol- 
low it, and a few residual samples are transient or drop 
to other fixed points mentioned below. As F becomes 
smaller, there appear three hierarchical branches (i) to 
the W state, (ii) to the 2-qubit Bell state (F = |) shared 
by two parties out of three, i.e., -^(|01) + |10))|0) Zl ' 2 ' 3 , 
where (Zi , Z2 , ^3) is a permutation of (A, B, C), or (iii) to 
the undistillable state x ( U P to local unitaries). Depend- 
ing on the initial entanglement, a branch is selected by 
the protocol. This hierarchy reflects the "onion-like" ge- 
ometry among different kinds of entanglement in 3-qubit 
mixed states |23j. As F approaches the "critical" region 
(F ~ 0.50 in Fig. |3J) for distillability, the characteris- 
tic number of steps toward every fixed point as well as 
the fluctuation of the progresses for different samples be- 
come larger. The fraction of the states which follow lower 
branches also increases. 

Conclusion.- Identifying a complementary (mutually 
unbiased) pair of stabilizer measurements, which replaces 
the conventional bilateral CNOT, as a key local oper- 
ation for distillation, we have proposed a 3-qubit W 
state distillation protocol. To our knowledge, it is the 
first protocol to distill directly multipartite non-stabilizer 
states. An extension to the n-qubit W state should 
be straightforward, introducing the general W basis by 
IfWbasis = 4= V™ Z 1 ■ ■ ■ Z^X 1 ! 1 * 1 Since our 

protocol distills a non-stabilizer state and stabilizer states 
on the same footing, our scheme may lead to a unified 
construction of direct distillation protocols for multipar- 



tite entanglement. It is still open whether a hashing pro- 
tocol can be made for non-stabilizer states without lo- 
cal depolarization which makes density matrices classical 
mixtures of pure states. Finally, since quantum comput- 
ers in which only stabilizer states are generated can be 
efficiently simulated by classical computers ^l)) the ap- 
pearance of non-stabilizer states, such as the W state, is 
necessary to exploit the power (universality) of quantum 
computers. Thus, the technique to purify such states, be- 
yond the "classical" parity check (exclusive OR via CNOt) 
for stabilizer states, might also give a new perspective on 
fault-tolerant quantum computation (cf. Ref. 24]). 
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